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Abstract 

The observed strong dark-to-luminous matter coupling ^ |51 13 suggests the existence of a some functional 
relation between visible and DM sources which leads to biased Einstein equations. We show that such a bias 
appears in the case when the topological structure of the actual Universe at very large distances does not match 
properly that of the Friedman space. We introduce a bias operator poM ~ Bp^is and show that the simple bias 
function h = l/(47rror^)6' (r — rmax) (the kernel of B) allows to account for all the variety of observed DM halos in 
astrophysical systems. In galaxies such a bias forms the cored DM distribution with the radius Rc ~ Ropt (which 
explains the recently observed strong correlation between Rc and Ropt [H), while for a point source it produces 
the logarithmic correction to the Newton's potential (which explains the observed flat rotation curves in spirals). 
Finally, we show that in the theory suggested the galaxy formation process leads to a specific variation with time 
of all interaction constants and, in particular, of the fine structure constant. 



1 Introduction 



The existence of Dark Matter (DM) has been long known 0]. It represents the most mysterious phenomenon of 
our Universe which stiU did not find a satisfactory explanation in modern physics. While more than 90% of matter 
of the Universe has a non-baryonic dark form, lab experiments show no evidence for the existence of such matter. 
The success of (Lambda) Cold Dark Matter (CDM) models in reproducing the large-scale structure is accompanied 
with a failure in describing the Universe on smaller scales. Indeed, it is now well established that in galaxies the 
Dark Matter (DM) density shows an inner core, i.e. a central constant density region (e.g., see Refs. [2] for spirals 
and Refs. for ellipticals and references therein). Such a feature is in clear conflict with A— CDM models which 
predict the presence of cusps (pdm ^ 1/?") in the inner regions of galaxies |5j (see however a more positive view in 
Ref. [S]). The situation is somewhat better for the Milgrom's algorithm [7] MOND (Modified Newtonian Dynamics). 
However, the existence of a very strong correlation between the core radius size Rq and the stellar exponential 
scale length Rd (or the optical radius Ropt^), Rc — 13 (5^^)^ °^ kpc, e.g., see Ref. P, rules out MOND as well. 
Indeed, according to Milgrom's algorithm the low acceleration regime triggers off at Rmond, when the gravitation 
acceleration g = GMgai/r^ drops below a fundamental acceleration oq 2 x 10~*cto/ (i.e., R%jond ~ GMgai/ao), 
and in general the two parameters Rd and Rmond are independent. By other words there should exist galaxies 
in which either Rd <^ Rmond, or Rd 3> Rmond- And indeed an example of such a galaxy has been recently 
presented in Ref. [H]. 

Thus we see that the modern theory of structure formation faces a rather difficult situation. Main alternatives 
to CDM, worm DM and self-interacting DM, seem to be ruled out by data on large scales (e.g., see Ref. [H] and 
references therein) , while the distribution of DM in galaxies rules out CDM ^ |2] and MOND . 

The correlation between the core size Rc and the optical size Ropt in galaxies of different morphological type ^ 
points clearly out to the presence of a very strong coupling between DM halos and baryons which surely requires some 
new physics. We recall that such a strong dark-to-luminous matter coupling (the so-called bias) is actually observed 
on all scales (e.g., Ref. [El EI)- In general this means the existence of a functional dependence or the so-called bias 
relation T^*^ = F^^{Tap) between DM T^/^ and the visible matter Tp^ sources. In the linear case the bias can be 
expressed by 

T^i^'^BT^,^ j B^, '^P{x,x')T^p{x')da', (1) 

J x' <.x 

where to save the causality the integration should be taken over the past-light-cone of the point x. In CDM 
models the bias relation appears as a result of the nonlinear dynamics during the structure formation and carries 
a nonlinear character, while on very large scales, where inhoniogeneities are still in the linear regime, such a bias 
should be viewed as the result of a generation process of primordial perturbations or merely as a result of the specific 
choice of initial conditions. In the present paper we consider the simplest case, i.e., the isotropic and homogeneous 
Universe with visible matter in the form of dust. Then the bias operator can be expressed via a single function 
B"^ [x, x') = {5p5^ + b{t,x — x'). Moreover, in such a case the bias function b{t,x — x') can be fixed from 

observational data, e.g., for Fourier transforms the bias relation gives 

T^^''{t,k)^B{t,k)T;mt,k) (2) 

which allows to find empirically the bias operator Bgmp- And it is quite obvious that the empirical bias operator 
Bemp (in virtue merely of its definition) will perfectly describe DM effects at very large scales (i.e. in the region of 
linear perturbations). Any actual specific source of DM (to fit observations) should reproduce properties of the bias 
operator B^mp in details. 

The bias relation allows to re-write the Einstein equations in the equivalent form 

- ^9^..R - 8nG {Tp, + Ff,,{T^f})) . (3) 

Now we can forget about the origin of the bias and study straightforwardly equations in the form Q . The advantage 
is that equations (jSJ do not imply the existence of any actual DM source. Therefore, with the same success we can 
interpret (jSJ as a specific modification of gravity. Most of modifications suggested (e.g., see Refs [ZIEUE]) can be 
reformulated in the form |j2Jl. In particular, for a point mass at rest equations © lead to a Modified Newton's law 

= + (4) 



^ Ropt is the radius encompassing 83% of the total luminosity of the galaxy. In the case of a (stellar) exponential thin disk Ropt is 3.2 
times the disk scale length Rd- 
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where in general the correction / {t, r) depends also on the position of the point source in space. We also note 
that such a modification can be equally interpreted as a specific "renormalization" of the gravitational constant 



In the present paper we discuss the bias relation which appears in the case when the topological structure of 
the physical space (i.e., of the Universe) does not match properly that of the Friedman space. It was demonstrated 
recently (e.g., see Refs. ^1^]) that in this case the standard Newton's law violates (there exist a range of scales 
ro < r < Tinax in which the gravitational potential has the logarithmic behavior, i.e., f {t,r) — r/rolnr). We 
show that the simple bias predicted in Refs. ^lEl b = l/Anro \r — r'\ 6 {r — r^nax) gives a rather good qualitative 
agreement with the observed picture of the Universe at smaller scales. In particular, such a bias allows to relate 
together a number of observational facts. Namely, the asymptotically flat rotation curves of spiral galaxies 
(which indicate that starting from some length scale rg the gravity force behaves as the cored distribution of 

DM density in galaxies the observed very strong correlation between Rc and Rd CO, and the fractal behavior 

in the distribution of galaxies (which has the dimension Z) ~ 2 and is observed at least up to 200Mpc QHI)- In the 
view of the modification of the Newton's law the last fact indicates that the maximal scale rmax after which the 
standard gravity law restores (e.g., it becomes F ~ again) should be rmax > 200Mpc [T7j). 

All these facts are well established and are beyond doubts. There were some debates in the literature about the 
fractal distribution of galaxies ^Hl- However, the test for the fractality is rather simple, e.g., if we consider any 
galaxy, surround it with a sphere of a radius R, and count for the number of galaxies N (i?) within the radius R, we 
find the law N (R) - R^ . And the value D is, in turn, not sensible to small perturbations of the galaxy distribution 
which may appear due to uncertainties in distances^. Moreover, the large-scale structure, e.g., the existence of huge 

100 — 200Mpc) voids with no galaxies inside and thin filled with galaxies walls (~ 1 — 5Mpc), is quite consistent 
with D ~ 2. Thus, it is safe to accept the fractal picture, at least up to 200 Mpc. 

2 Origin of the bias 

In the present section we show that a non-trivial topological structure of the physical space can quite naturally give 
rise to the origin of the bias |13l I14| . Indeed, in considering astrophysical systems we use an extrapolation of spatial 
relationships which are well-tested on considerably smaller scales. Therefore, if the topological structure of the actual 
Universe at very large distances does not match properly that of the Friedman space (the open, fiat, or closed model) 
we naturally should observe some discrepancy. To describe such a discrepancy we first consider an example from 
solid state physics. 

Consider a medium of a low density at very small temperatures. From the thermodynamics we know that most 
of systems at a sufficiently small temperature acquire a crystal structure. However, in actual systems such a crystal 
has never an ideal character but includes different distortions. Moreover, when a system has a rather low density 
and the rate of freezing is rapid enough, such a system will include considerable voids and the spatial distribution 
of particles in the system acquires, in turn, quite irregular character. Elementary excitations (or quasiparticles, e.g., 
electrons of the conductivity, phonons, etc.) in the given system do exist only within the crystal and from their point 
of view the physical space (the crystal) possesses a rather non-trivial topological structure. From the mathematical 
standpoint the non-trivial topological structure can be accounted for as follows. 

Consider a volume V in R'^ occupied with a system and let H be the Hilbert space for a free particle (the space 
of functions on V). Let {gk (x)} { x £ V) be an arbitrary basis in H. Physically, the basis represents a set of 
eigenvectors for a complete set of observables. E.g., for a scalar (without the spin) particle we can use the coordinate 
representation (i.e., gkix) = 6 {xk — x) is the set of eigenvectors for the position operator Xg^ = Xkgk, Xk € V) 
or the momentum representation {gk{x) — {V)^^^"^ exp (ikx), so that Pgk = kgk). The basis is supposed to be 
normalized {gk,gp) — ^kp and complete 5fe (2^) 5fc (^') = 5{x — x'), where x, x' G V. The fact that our system 
has an irregular distribution in V (i.e., V includes also voids) means that some states in H cannot be physically 
realized for particles of the system (at least for small temperatures when the structure of the crystal does not change) . 
Thus, we have to restrict the space of states H to the space of physically admissible states Hphys = KH, where 
K = {KY is a projection operator. In the basis of eigenvectors the projection operator K takes the diagonal form 
{fi,Kfk) — Kik — NkSik with eigenvalues Nk — 0, 1. Thus, an arbitrary (but physically realizable) state of a particle 
is biased and can be presented as ^phys = K-^-^^ip = ^ ^/NZakfk {x). Thus we see that topological structure of the 
system is described by the bias (projection) operator K. In particular, all physical observables acquire the structure 

■^The misunderstanding may appear if one performs an averaging over the central position of the sphere in space. In this case one gets 
nothing but the trivial result D 3. 
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Ophys = K^/^ O K^/^, while the physical space Vphys of the system represents the space of eigenvalues Xk € Vphys 
of the biased position operator of a particle Xphys = K^^^ X K^/^. 

In the example described the bias operator is diagonal in the coordinate representation (i.e., Nk = 0, when Xk 
belongs to voids and = 1 as Xfc belongs to the crystal). However, we can also consider a more general case when 
K and X do not have common eigenvectors (i.e., [K, X] 7^ 0). In the last case the spatial structm'e of the crystal 
remains unspecified. This means that in such a system the position operator cannot be a good observable (at least 
while the topological structure of the system conserves, i.e., Kik = const, which is always fulfilled at sufficiently 
small temperatures). We also note that from the point of view of the mathematical coordinate space (i.e., R^) 
the space Hphys is not complete, i.e., J2 ^kfk i^) fk {x') = K{x,x') = K^f^S {x - x') K^/"^ ^ 5{x-x'). Thus, 
we see that the function K{x,x') plays here the role of the delta function. And only in the case when both K 
and X can be diagonalized simultaneously the biased delta function K{x,x') reduces to the ordinary delta function 
K{x, x') = 5 (x — x') 9{x, Vphys), where 9{x, Vphys) is a characteristic function, i.e., 9{x, Vphys) = as x ^ Vphys and 

0{X, Vphys) = las X € Vphys- 

At very low temperatures the structure of the crystal conserves. This means that the projection operator K 
represents an integral of motion (commutes with the Hamiltonian of the system). Therefore, we can state that 
elementary excitations (quasi-particles) represent eigenvectors for the projection operator i.e., the wave function of 
an excitation can be expanded as 'ipphys — V^kO-klk (x), while the energy of the system can be represented as 

E^J2^keik)a+ak, (5) 

where s (fc) is the energy of a quasi-particle. Thus, we see that the non-trivial topological structure of the system 
defines the measure (i.e., the density of degrees of freedom) which can be accounted for by the formal substitution 

k k 

(indeed, the algebra of physical observables modifies as A = BC Aphys — BphysCphys — K^^^ B KC K^^^ 
and {B KC)ij = -^kBikCkj)- Any point source for quasiparticles is always biased (as compared to the simple 
topology case), i.e., acquires a specific distribution in 

S{x^ x') ^ K{x, x') = K^/^ {x - x') (7) 

which reflects the topological structure of the system (the discrepancy between Vphys and V). In particular, the 
actual physical volume occupied by the crystal is given by Vphys = Jy ^ d^xcfix' ^ V . 

The above construction generalizes straightforwardly onto relativistic particles. In a curved space the one-particle 
Hilbert space is not well defined, for particles are actually not free. This means that in general there is no such an 
observable as the position operator X or the momentum P to classify quantum states. We recall the well-known 
fact that even in the flat space the momentum of a particle can be considered as a good operator, while the position 
operator is not. It can be defined though by means of the Newton- Wigner construction 19 . Thus, in this case the 
space of quantum states is constructed as follows. 

Consider an arbitrary set of solution to the wave equation'^ 

(o + ^R + TO^^ fk = 0, (8) 
(where Dfk — y=9a {\/~~gg°'^dpfk)) which obey the normalization conditions 

{fk, /,) - - (A*, /;) - 4„ (At, /,) = , (9) 

and the scalar product is defined as (e.g., see Ref. |2(J|) 

(/i, /2) = ^ j (/r {x) [x) - h (x) V^/r (x)) V^dS^ (10) 

Then the space of one-particle quantum states is defined as the space of "positive frequency" solutions {fk}- 
And again in simple cases a non-trivial structure of the physical space can be accounted for by the fact that some of 
one-particle quantum states cannot be physically reahzed, i.e., we should project the space of states to the space 

^If we require that the topological structure should be invariant under conformal transformations, then we should set m = in ISl . 
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of physically admissible states H^f^^^ — KH^ . In general the projection (bias) operator distinguishes a particular 
(preferred) basis {fk} in terms of which it can be presented as^ 

(x, x') - 5] iV, (/, {x) ft (^') - ft {x) fk ix')) , (11) 

with eigenvalues = 0, 1. Thus, physical fields can be defined as biased fields 

0ph,. = kI'^<P - 5] (flfc/fc {x) + 4fl {x)) (12) 

and the nontrivial topological structure of space is reflected in the fact that some modes never enter the expansion 
ltT3 (i.e., for which =0). And again any physical observable (i.e., every operator) can be expressed as <Dphys — 
K^/^ O K^/^. E.g., in the case of a scalar field the mean value for the stress energy tensor is biased as 



rpphys 
^af3 



nk)^{nk k^'^T^pK^'^ Uk) = ^ TVfe (1 + 2nk) T^p [fk {x) , fl {x)\ , (13) 

A; 

where Tap [0, 0] is given by the bilinear form 

Tap [0, 01 = (t>,a(t>:p - ^5a/3 - m^(t>^*) (14) 

and \nk) = 11 ("fe')"^^^ (^fc)"*" 1^^)- Green functions for the physical scalar field (e.g., Feinman propagator 

iGp {x, x') — (0 \T(j)phys (x) (f>phys ix')\ 0)) obey formally the standard equation 

(□ + m^) Gf (x, x') = A {x, x') ^ - {~gy^'^ 5{x- x') . (15) 

However the r.h.s. of this equation is not the delta function any more but physical or biased delta function Q (i.e., in 
terms of the coordinate space it acquires an additional distribution in space A (x, x') = —K^^'^ {—g (x)) S {x — x') 

{—g {x'))"'^^'^ K^^"^). In this manner we see again that the role of the bias operator (and that of the structure of the 
physical space) is the specification of the density of degrees of freedom 10. 

In conclusion of this section we point out to the two important facts. The first is that the bias (|ll|l includes a 
non-linear dependence on the metric gap via the solution of Eq. ||SJ|. And the second is that the projection operator 
(bias) discussed above restricts strongly the topological structure of the physical space. Indeed, by the construction 
the projection K — (K)'^ means that the physical space Vphys represents a subspace in (i.e., Vphys C or 
in cosmology it should represent a subspace of the Friedman space). In the most general case however such an 
embedding may not exist. By other words an arbitrary physical space (of an arbitrary topological structure) cannot 
be projected to the Friedman space (or R^) without self-intersections (i.e., K ^ {K}'^) ■ This, in turn, leads to a 
generalization of the bias operator Hll() to the more general case (e.g., see Refs. ^J^J) which naturally leads to 
the generalized statistics of particles. From the formal standpoint such a generalization is expressed by the fact that 
eigenvalues Nk of the bias operator K can be arbitrary integer numbers Nk = 0, 1, 2, (A^l 7^ Nk and K 7^ (K)'^). 

To illustrate the last statement we can consider an example from solid state physics. Suppose that the system 
discussed in the beginning of this section has locally a two dimensional character (i.e., locally Vphys represents a 
two-dimensional crystal). Then we can attempt to describe such a system in terms of R^. If we project Vphys onto 
R^, then we find that in the case of an arbitrary topology of the two-dimensional crystal Vphys the bias operator 
will have eigenvalues Nk = 0,1,2,.... E.g., if K is the diagonal in the position representation (i.e., [KX] = 0), 
then Nk is merely the number of different points of the crystal (i.e., the number of two-dimensional sheets) which 
correspond to the point Xk G R^- All such points however have different positions in R^, i.e., they differ in the extra 
coordinate z^{a = 1,2, ...,Nk) orthogonal to R^. However, if the Hamiltonian of the system does not include the 
extra coordinate Zk, it is not measurable (without additional means) and states, which differ in the extra coordinate 
only, become physically indistinguishable and quasi-particles will obey a generalized statistics. In particular in the 
given example Nk gives the maximal number of electrons which can occupy the same position Xk G R^. For more 
details see also Refs. [THI^ . 

In this manner we see that a non-trivial topological structure of the physical space (as compared to the coordinate 
space) can indeed produce a specific bias of all observables. We note that in this case the field theory does not change 
at all, i.e., the mathematical structure of equations of the motion (e.g., the Einstein equations) remains the same. 

■^Wc note that the operator defined in llli acts in the space of fields if) (x). In the one-particle Hilbert space it has the standard 
form i? = EA^fc life) (Ifel- 
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What is actually modified here is spatial properties of physical fields^ which are expressed by expansions of the type 
l)12|l . In particular, every discrete point source (e.g., a galaxy or a star) is not the Dirac delta function any more but 
acquires a specific distribution in space (e.g., see (T))) which reflects the topological structure of the physical space 
(the density of degrees of freedom Nk). 



3 The bias function b (r) 

In what follows we, for the sake of simplicity, restrict our consideration to the Newtonian limit (for the range of 
applicability of this limit see, e.g., Ref. 22 ). In a homogeneous and isotropic Universe the set of solution @ can 
be taken in the form fk = {2n)^^^^ gk {t) e^^"^ (i.e., states of particles can be classified by wave numbers fc), while 
the density of states is an arbitrary function of |fc|. If we assume that topology transformations have stopped 
after the quantum period in the evolution of the Universe, then the function will depend on time via only the 
cosmological shift of scales, i.e., kit) ^ l/a{t) (where a{t) is the scale factor). Thus, any point source undergoes 
the bias 

oo 

2^ / y^'^'' ) -^^T 



6{r) ^ A(f,t) = / {N,k') (16) 





The case of a simple topology corresponds to Ni^ — 1, while in a non-trivial case {Nk — 1 7^ 0) every point mass 
Mo is surrounded with an additional spherical " dark" halo 

00 

/ N , ^ , / X Afo /" , , , , , , ^sin ikr) dk „, 
RDM (r,t) = Moh{r,t) = ^ j [Nj, (t) - l)k^^-l— (17) 



and the Newton's potential modifies as 

<t> = ^(l + /(r,t)), (18) 

where the correction / (r, t) relates to the bias function b (r, t) according to (/ (r, t)' — df /dr) 

h{T)=-l^. (19) 

Thus, the relation between visible matter pms and DM is indeed given by 1^ which in the Newtonian limit for the 
homogeneous and isotropic Universe reduces to 

PDM (r, t) = Bp^^s = j b{\r-f'\,t) p^^s {r, t) dV . (20) 

The explicit specification of the bias function b (r, t) is, in the first place, the problem of observational cosmology. 
Indeed, for Fourier transforms there is a linear relation between DM and visible sources 



Pdm 



(fc, = b {k, p,,, {k, (21) 



which allows to find empirically the bias operator Bemp (we recall that the total source ptot ~ Pdm + Puis can be 
restored from the measured spectrum of AT/T in CMB and the observed peculiar velocity field). It is quite 
obvious that such an empirical bias operator Bemp (in virtue merely of its definition) describes perfectly DM effects 
at very large scales (where inhomogeneities have the linear character) . The nontrivial moment here is that all theories 
which predict the same bias b (r, t) for the modern Universe are observationally indistinguishable (at least it requires 
involving more subtle effects). We also note that in the more general case the bias relations should be described 
by two functions pdm — bpPms and pdm = &pPuis (where p is the pressure) which for a homogeneous distribution 
reduce merely to functions of time b'^ ^ (t). Thus, the bias relation give the possibility to account for Dark Energy as 
well (i.e., the observed^ accelerated expansion of the Universe 

A specific feature of CDM models is that the relation between the two sources appears as a result of the dynamics 
and, therefore, the effective bias function b (r, t) carries in general a nonlinear character. The "great" success of CDM 

^In general physical fields should be understood as generalized fields Refs. BH - 

®We point out however that the accelerated expansion can not be considered as an established fact yet, for the presence of considerable 
uncertainties of a theoretical character. 
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models in reproducing the large-scale structure (LSS) of the Universe is somewhat exaggerated, for at very large scales 
density perturbations are still at the linear stage of the development and, therefore, the bias fcemp (j^, straightfor- 
wardly defines the set of appropriate initial conditions bemp (k, — D (t) bo (jtj (where &o (k) — p^j^.j ^fc^ /Pvis (k^ 

and D (t) accounts for the evolution of perturbations) depending on the exact behavior of the scale factor a (t) . In 
this sense LSS alone in principle cannot distinguish a model. On the contrary, at smaller scales (e.g., in galaxies and 
clusters) perturbations are in a strongly nonlinear regime, the bias operator B acquires a nonlinear dependence on 
the distribution of matter and CDM models fail 2 , 3 . 

Leaving the problem of the empirical determining of B aside, in what follows we consider a model expression for 
the bias b (r) 



b (^) = 7ZrT2 ^ - ) ' (22) 



1 

where 6 (x) is the step function, b (r) produces the correction to the Newton's potential H18|l of the form 



f f \ _ I r/ro\n (r,„axe/r) , as r < 

' max; /oo^ 

' max /ro, as r > 

' max • 



Such a bias was derived in Refs. |13II14) for the case of a homogeneous and isotropic Universe under the assumption 
that the topological structure (i.e., the number density of degrees of freedom Nk) of the early Universe is described 
merely by the thermal equilibrium state^. Presumably, topology changes have occurred during the quantum stage of 
the evolution of the Universe and at present are strongly suppressed. This means that after the quantum period the 
topological structure remains constant. Therefore, the isotropic cosmological expansion is accompanied only with 
the cosmological shift of the parameters ro and rmax (i-e., ro^max (^) = a (i) Hj.max) without any change in the form 
of the bias function (|22|l . 

After the radiation dominated stage, however, the small initial adiabatic perturbations (which are directly mea- 
sured in CMB e.g., by WMAP '^) start to grow and considerably shrink the Universe from galactic to supercluster 
scales. The latter results in the further transformation of the bias function b{\x — x'\) ^ b{\x — x'\ , x' , t). To derive 
rigorously the bias in a general inhomogeneous case we have to construct a set of exact solutions to the wave equation 
^ which in turn depend on the distribution of matter and, therefore, on the bias. In the simplest case however 
the inhomogeneity of the Universe can be accounted for by an additional dependence of the parameters of the bias 
function H22|l on the position in space. Indeed, the adiabatic growth of density perturbations can be viewed as if the 
rate of the expansion were different in different parts of the Universe a{t) a {t, x) which produces the respective 
shifts ro,max(i,a;) ~ a (i, a;) ro^max- Such an additional shift is considerable indeed, e.g., the mean density of our 
Galaxy has the order pg ^ 10^ Pcr (while the density behaves as p ~ l/a"^) and therefore for our Galaxy r^o should 
be less in 10^ times, than the respective mean parameter ro for the homogeneous Universe. 



4 The bias function and Dark Matter halos 

It is rather surprising that already the simplest function (|22|l shows a rather good qualitative agreement with the 
observed picture of the present Universe. First of all it is consistent with the observed cored distribution of DM in 
galaxies Indeed, if pms (r) is a rather smooth monotonously decreasing function of r, then from H22|l and (|20() 

we find that DM density reaches the maximal value in the central region of a galaxy (i.e., as r ^ Rjj, where Rd has 
the order of the stellar exponential scale length) 

PDA! (r) =i PDM (0) = / J^-^Pv^s (^, t) dV , (24) 

J 47rror'^ 

while for r 3> Rd we find poM {r) — M^s/ (47rror'^) (where Myis = / pmsdV) which can be combined by the 
interpolation formula 

R^ 

PDM {r) = PDM (0) 2^2 ' (25) 

+ r 

where Rq = Myis/ (47rro/ODj\/ (0)) — a^Rj^) ) which explains the observed strong correlation between Rc and Rd 
We note that the actual value of the ratio Rc / Rd — ol depends on the distribution of the visible matter in a galaxy 
Pvis {fjt) and the definition of Rd (e.g., if we assume in (|^ that pyis = p within the baU r < Rd, then = 1/3). 

'^We note that the actual bias depends on the specific picture of topology transformations in the early Universe and may differ from 
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The bias (|22|l shows also that in the interval of scales r < Tmax the dynamical mass of a point source increases with 
the radius as Mdyn = Mq (1 + r/rp), while for r > rmax it acquires a new constant value Mmax ~ Mq (1 + rmax/'^o) 
and the ratio Tmax/fo defines the fraction of DM in the total (baryons plus dark matter) density. 

The minimal scale tq is different for different galaxies (i.e., = rg (x) is a slow function of the position) and it 
has the order ~ 1 — 5 kpc (it is the scale at which DM starts to show up), while the value of Tmax is not so well 
fixed by observations. The analysis of the mass-to-light ratio M/L shows that it increases with scales for galaxies 
and groups but flattens eventually and remains approximately constant for clusters (e.g., see Ref. |25)). This gives 
an estimate rmax ^ 1 — 5Mpc or Tmax/^'o ^ lO'^- Such a fraction of DM is indeed observed in LSB (Low Surface 
Brightness) galaxies in which the ratio can reach M/L ^ (200— 6OO)Mo/L0. It however looks inconsistent with 
predictions of CDM models and observed peculiar velocities in clusters which favor pdm / Pb ~ 20. The most drastic 
estimate comes from the observed fractal distribution of galaxies which suggests r,„ax ^ 200Mpc and rmax/^o ^ 10^ 
[T7|. We however note that the absolute boundary for r^ax is given by the Hubble radius Tmax < Rh which gives 
Tmax/ro < Rn/fo ~ 10^ — 10'', while all values Tmax > Rh are indistinguishable from observations. 

It turns out however that all those estimates are consistent with each other and give only the lowest boundary 
for the DM fraction, for in any system some essential portion of DM forms an inner core (i.e. the central constant 
density region) and does not contribute to the local dynamics. Indeed, DM consists of spherical halos (|17() around 
every point source and, therefore, the relationship between the baryon density and DM has a non-local nature with 
the characteristic scale Tmax- The density of DM in a point of space (and respectively the local dynamics) is formed 
by all sources within the sphere of the radius Tmax and it depends essentially on the distribution of the sources. E.g., 
if we take pms (x, t) — Ma5 {Ra), then from and we get for DM density 

/ ,\ ^'-^a ^ ''max {Pvis) /r,c:\ 

where Ra = \x — Xa {t)\ and {pvis) — Sfl„<r„„x ^^"Z (4/37rrjJjax) is the mean density of the visible matter within the 
sphere of the radius rmax. For a uniform distribution of matter this reads (pdm) — (j'max/j'o) (pvis)- From H26|l we see 
that the DM density reaches the minimal possible value 1/3 (pdm) in the case when all sources are at the distance 
Ra = ''max (e.g., in the center point of a void), while according to H25|l at a source Ma it has a local maximum 
Pdm — /'"o) {Pa) /3 (where £a is a characteristic dimension of the source and {pa) = 3Ma/47r€^ ). 

(|26|) shows that DM halos smooth the observed strong inhomogeneity in the distribution of baryons which con- 
siderably reduces the inhomogeneity in the total density. By other words a considerable portion of DM acquires the 
cored (|25f) (i.e., the quasi- homogeneous) character and switches off from the local dynamics. This, in turn, leads to 
a renormalization of the maximal scale r,„ax ^ R* in (|22|) and, therefore, changes the fraction of DM observed in a 
system p'^j^j/ Pb ^ R*/to- In such a picture the scale is a specific parameter of a system and this explains the 
small value for the ratio i?*/ro observed in clusters. 

Indeed, consider a group of galaxies of the characteristic dimension L. Such a group can be characterized by 
the mean density {pdm) = ^/L^ J Pdm {x,t)(Px and perturbations Sdm {x,t) = pdm {x,t) / {pdm) ^ ~ 1- Near a 
particular galaxy in the group (r^ {t) = and Mg <C X) ^^a) we find from 

5dm {r)^^~ 1, (27) 



where i?* is the effective size i?* of the DM halo 



Rl _ Mg 



47rr3 {pdm)i 



(28) 



For r > we see that Sdm < and in the interval L > r > R,, this function oscillates around the zero point (the 
exact behavior depends on the distribution of galaxies in the group and is not important). 

The homogeneous background contributes only to the local Hubble flow which can be accounted for by the 
expanding reference frame x = a{t)r (e.g., see Ref. (22,). Thus, the actual Newton's potential of the galaxy takes 
the form 

(50 (r, t) = -Go" + SFdm (r, t)) , (29) 



with 

SFdm (r, t) - ^ m / '^ijl?'^ + 2/3n {pdm) = Mgl^^^ + p (r, t) (30) 
\r r.j [t)\ r 
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where we subtracted the homogeneous component S(f> = (fi— {(f>) (with (0)^ = 2/3TTGa^ {p)^ r^) [21], (r, t) accounts 
for variation of 6dm for r > R^, and / (r, i?,) is given by (|23|l with the replacement rmax ^ The function SFdm 
defines the contribution of the DM halo and we recall that the use of the empirical bias function b^mp (j", '''ji) (or 
equivalently f {r.t)) automatically reproduces all actual DM halo s in astrophysical systems. 

Thus, we see that near a source the function 5Fom has the logarithmic behavior^. At the distance i?, the 
logarithm switches off and the ratio R*/ro defines the maximal value for the DM mass in a galaxy or a cluster 
which can be observed from the local dynamics. We recall that the value rp is different for different galaxies. In 
addition to this fact the expression 1)28(1 shows the general tendency that the ratio R*/ro (and therefore the maximal 
discrepancy between the dynamical mass and luminous matter) is smaller in high density regions of space and larger 
in low density regions. This qualitative feature agrees with discrepancies observed in LSB and HSB galaxies. 

5 The background distribution of baryons and rmax 

Consider now properties of the homogeneous and isotropic background. In the standard models there exist the only 
case which corresponds to the homogeneous distribution of baryons. If we accept the bias of wave equations H15|l. 
there appears a new possibility. Indeed, the homogeneity of the Universe (or the cosmological principle) requires the 
total distribution of matter (baryons plus dark halos) to be homogeneous, while properties of the baryon distribution 
are not fixed well. The latter may have a quite irregular character. Exactly such a situation takes place in the case 
of a fractal distribution of baryons. Consider a sphere of a radius r. Then the total mass within the radius r is given 

by 

Mtot (r) ~ mfc (1 + r/ro) m (r) + SM (r) , (31) 

where Nb (r) is the actual number of baryons, nib is the baryon mass, and 6M (r) accounts for corrections and, in 
particular, for the contribution of dark halos of baryons from the outer region. The homogeneous distribution means 
that the total mass behaves as Mtot {r) = (p) V (r) ^ r^ . And for r 3> tq this can be reached by the fractal law 
N (r) ~ r^ with D k, 2 (the exact equality cannot be reached, for the presence of the additional term 5M (r)). Such 
a law works up to the scale Tmax upon which the distribution of baryons crosses over to homogeneity^. 

There exists at least two strong arguments in favor of the fractal distribution of baryons. The first argument is 
that the fractal distribution is more stable gravitationally. Indeed, let us fix the total density O.tot = Ptot/ Per ^ 1 
(where pcr is the critical density) and the baryon fraction pb/ Ptot ~ fo/fraay.- In the case of the fractal distribution 
this fraction reaches only at scales r > r^ax, while at smaller scales baryons are distributed rather irregularly. 

Consider first a homogeneous distribution of baryons. Now if we consider a small displacement of a particular 
baryon (or of a homogeneous group of baryons), then such a displacement will produce the same displacement of 
the dark halo (attached to the baryon). So the resulting perturbation increases in rmax/'^o times. The maximal 
scale r„iax should be larger than 100 — 200Mpc, and therefore the increase should be more than 10^5 - 10'5. In the 
primordial plasma the domination of radiation prevent the growth of perturbations in the gravitational potential 
and, therefore, such fluctuations are strongly suppressed. However, there also do exist collective fluctuations in the 
density of baryons which do not affect the metric perturbations^^ and the total density of matter. According to 
lfTK)l such fluctuations do not affect the total (effective) charge density and, therefore, the radiation dominated stage 
cannot prevent a specific redistribution of baryons. By other words perturbations of such a type could increase 
long before the recombination. They do not change the total density Sptot = 5pb + SpoM — const and can be 
called compensational sound waves. In the very early Universe high temperatures transform baryons from the more 
constrained state which corresponds to a homogeneous distribution of baryons to the less constrained and more stable 
state which corresponds to the fractal distribution. We note, however, that during the radiation dominated stage 
when 6pb + SpoM ~ and Qb + ^ dm ^ 1 perturbations in the baryon number density cannot grow to an arbitrary 
large value, but are restricted by Sib < 1, (i.e., 6pb/pb < Ptot/pb ^ fmax/''o)- 

Consider now the case of the fractal distribution. According to H31|) the fractal distribution of dark halos forms 
the homogeneous background of the total density. Now any small displacement of a baryon does not change the 
character of the dark halos distribution and, therefore, the increase is essentially suppressed (rmax/''o ~^ R*/fo)- By 
other words the stable equilibrium distribution can be defined as such a distribution of baryons for which perturbations 
in the baryon density produce the minimal response in the total density. The bias of the electromagnetic field ((15(1 

^We note that in the presence of a continuous medium (e.g., of gas) the behavior may essentially change. 

^The distribution of stars in galaxies shows also a fractal behavior. In this sense we can say that the fractal law forms the cored 
distribution I25i with ~ Rh- 

^"The presence of metric perturbations at some level Aptot/ptot ~ 10^^ is essential however, otherwise the fractal structure in baryonic 
matter will not form. For fluctuations the bias relation reads ApoM = BAp^ig — {BAp^ig) = FAp^is, which defines a new operator F. 
Thus, fluctuations which obey (_F — 1) Ap„is fa do not affect the metric and Aptot ~ const. 
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insures the absence of strong fluctuations in the CMB temperature caused by the fractal distribution of baryons. 
This may be used to estimate the vahie of the fraction r„iax/?'o- 

Indeed, the first estimate comes from the upper boundary for the scale of the cross-over to the homogeneity 
in the observed galaxy distribution Tmax > 100 — 200Mpc which gives rmax/^'o > 10^. From the other side, the 
observed CMB gives AT/T — ^Aptot/ Ptot ^ 10~^ at the moment of recombination, and the fractal distribution 
causes perturbations in the total density Aptot ^ {R*/ro)pb (where the factor i?*/ro appears as the contribution 
from dark halos) and therefore Aptot/ptot > {R*/ro) Pb/ Ptot ~ R*/rmax < 10~^. We see that both estimates agree 
and give Tmax/'^o ^ (^*/'"o) 10^- As it was shown above the ratio i?*/ro takes the minimal value for the equilibrium 
fractal distribution. So that the value fmax (which is the scale of the cross-over to the homogeneity in the visible 
matter) will increase, if at the moment of the recombination the ideal fractal distribution had not been achieved yet. 

The second argument is based on a more correct interpretation of the dark matter effects. Indeed, the bias 
of the wave equation fl^ should be understood as the fact that at large scales our Universe possesses a rather 
unusual geometric (or topological) properties. These geometric properties are reflected in the behavior of the Green 
function l|15fl which for r > tq acquires effectively the two-dimensional character (e.g., for Nk ~ 1/ {kro) we get 
G (r, t) ~ ^ In {t — r) / {t + r)) and, therefore, such a geometry should be reflected in the distribution of matter 
(sources). By other words at scales r > ro our Universe acquires an effective dimension D « 2 (e.g., see Ref. |14| ') 
which explains the two-dimensional character of the spatial distribution of baryons. By other words we may imagine 
that our Universe represents a fractal (the space is "more dense" on a fractal set than outside (e.g., see Ref. 52)) 
and within such a fractal the matter has a homogeneous distribution. In such a picture the fractal distribution is the 
only thermal equilibrium state. We note that in the case r^ax < oo such a state can never be utterly homogeneous 
but always includes equilibrium fluctuations of the order Aptot/ Ptot ^ fo / fmax- 



6 Variation of interaction constants 

In the present section we show that the structure formation in the present Universe leads to a specific variation with 
time of all interaction constants. As an example we consider the variation of the gravitational constant. Indeed, the 
cosmological evolution is described by the scale factor a (t) which obey the equation (22j (we consider the case p = 0) 

(Pa AttG 47rG' / r^axN , > .„„^ 
— :^{Ptot)a = 5—1 + ]{pvis)a. (32) 



This equation can be interpreted as if the gravitational constant renormalizes as G — > G = G (1 + r,nax/?'o) (we recall 
that in the inhomogeneous case it depends on scales as well). 

The mean density of the visible matter behaves as (pvis) ^ l/a^. Thus, the evolution of the scale factor a{t) 
depends essentially on the behavior of the ratio rmax/'ro- During the radiation dominated stage (ptot) = (p-y), the 
growth of density perturbations is suppressed and therefore the bias function H22|l in the commoving frame (i.e., in 
the expanding reference system x = ar) does not change. Thus, during the radiation dominated stage the ratio 
Traay^/fo — const. This remains true and on the subsequent stage, while inhomogeneities in the total density remain 
small d — Aptot / Ptot ^ 1- The situation changes drastically when the inhomogeneities reach the value ^ ~ 1. Upon 
this moment the time shifts of the two scales vq and rmax disagree. Small scale inhomogeneities develop first and 
switch off from the Hubble expansion. This leads to the monotonic increase of the effective gravitational constant 
G, i.e., of the ratio rmax/''o ^ a^, which gives for the matter density {ptot) ^ a^~^. While inhomogeneities remain 
small S < 1, both scales increase with time as ro, ^max ~ a, and the exponent /3 ~ 0. When S reaches the value S > 1 
the scale vq starts to collapse (galaxies start to form), while rmax is still increasing rmax ~ 0,. This leads to the fact 
that the exponent becomes /3 > 1 and DM behaves as "Dark Energy", e.g., in the case /3 = 3 DM behaves as the 
negative Lambda term A = —AttG (pdm) = const. This kind of regime ends either when the collapse of the scale ro 
ends (galaxies have stabilized and ro ~ const and G ~ a), or when the maximal scale rmax sufhciently deviates from 
the Hubble law r^ax ~ a- 

1/3 

The behavior of the minimal scale ro follows the local dynamics and can be estimated as ro '--^ arg, where 
Sq is the mean perturbation within the radius vq and the parameter ?o = const. Analogously, the maximal scale 
is given by rmax ^ i5max'^ar,nax, which gives rmax/^^o ~ {Sma.x (t) /Sq (t)) ^^'^ and therefore the effective gravitational 
constant depends on time as G {t) ^ G (^l + ((5max (t) /So (i))^^^'^^ • 

The fact that the bias operator reflects the topological structure of space means that all interaction con- 
stants undergo an additional renormalization (e.g., see Ref. j23) and acquire the same dependence on time. 
E.g., the fine structure constant takes the form a{k,t) = b{k,t)a which gives for homogeneous fields a{t) « 
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a yl + ((5max {t) /6o (t)) ^^^j. It is remarkable that a small variation of the fine structure constant seems to be 

observed at high red shifts |26| . 

In conclusion of this section we note that the decrease of the scale tq (t) during the structure formation can also 
be used to explain the apparent acceleration of the Universe which seems to be required by observations of the type 
la supernovae ,24, ■ Indeed, according to ((T3|l . lfTB|) . and at large distances r ^ tq the Green functions behave 
as G ~ l/ro and therefore the apparent luminosity will also behave as L ^ Lq/tq. Thus, the decrease of the scale 
To will formally look as a very strong evolutionary effect E = L/ L ^ —'f'o/fo > , which produces correction q — s- 
q'^ff = q — E/H, e.g., see Ref. |57] (where H = a/a and q = —(Pa/dt^/ {aH^^). Thus, the observed acceleration 
g < may merely mean nothing but the strong evolutionary effect caused by the variation of rg. 



7 Conclusion 

In conclusion we briefly repeat basic results. First of all from the observed strong dark-to-luminous matter coupling 
PmO] we derive the existence of a bias relation T^^^ = F^^^ which allows us to re- write the Einstein equations 

in the equivalent biased form R^j.v — \g^ji.vR = SttG {T^^ + F^i,{Taf})). The biased Einstein equations straightforwardly 
predict the presence of a specific correction to the Newton's potential for a point source = ~GM {1/r + F{r, t)). 

The bias may have an arbitrary nature, CDM, MOND, any modification of gravity, etc., which does not change 
the phenomenological results of this paper. We however have suggested the bias which naturally appears in the case 
when the topological structure of the actual Universe at very large distances does not match properly that of the 
Friedman space (the open, flat, or closed model). In that case not only the gravitational potential but also all other 
physical fields undergo the bias and display some discrepancy (i.e., the presence of DM halos around every point 
source S {x — x') A{x — x')). 

In the linear approximation the bias relation pdm — Bp^is is described by the function h{r,r' ,t) (the kernel of 
the bias operator) which admits the empirical definition. Then 6emp (r, r', t) (or equivalently its spectral components 

) gives a rather simple tool for confronting a theory of the structure formation with observations. Any 

acceptable theory has to reproduce in details the specific form of the bias function h^mp- 

We have demonstrated that a specific choice of the bias 6=1/ (47rror^) 9 {r — r^ax) (which is predicted by 
topology changes in the early Universe shows quite a good agreement with the observed picture of the modern 

Universe (e.g., the fractal distribution of galaxies, cored DM distribution in galaxies and rich clusters, variety of DM 
halos, etc.). It however considerably changes the estimate for the mean density of baryons {pdm) / {pvis) ^ ''max/''o 
(this in turn is not in a conflict with observations, for in the standard models the estimate Vth ^ 0.05 is model 
dependent and uses essentially the idea of the homogeneous distribution of baryons). 

Finally, we have shown that the galaxy formation process causes a decrease of the minimal scale rg {t) (and the 
increase of the ratio ■Tmax/^'o) and this gives rise to a specific dependence on time for all interaction constants. In 
particular, this may give an explanation to the observed variation (a small increase) in the fine structure constant 

M- 

In conclusion I would like to acknowledge discussions with D. Turacv and M. Milgrom on problems of modified 
gravity which inspired me for this research. I indebted to P. Salucci for attracting my attention to the recent 
observational results in DM physics 000]. I would also like to acknowledge the hospitality of D. Turaev and the 
Center for Advanced Studies in the Ben Gurion University of the Negev where this research has been started. 
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